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ABSTRACT 

We revisit the concept of a blazar sequence that relates the synchrotron peak frequency (i^poak) in 
blazars with synchrotron peak luminosity (Lpoak, in i^L^) using a large sample of radio-loud AGN. 
We present observational evidence that the blazar sequence is formed from two populations in the 
synchrotron ^'peak — Lpeak plane, each forming an upper edge to an envelope of progressively misaligned 
blazars, and connecting to an adjacent group of radio galaxies having jets viewed at much larger 
angles to the line of sight. When binned by jet kinetic power (Lkin! as measured through a scaling 
relationship with extended radio power) , we find that radio core dominance decreases with decreasing 
synchrotron Lpeak, revealing that sources in the envelope are generally more misaligned. We find 



population-based evidence of velocity gradients in jets at low kinetic powers (^ 10 



42 



10 



44.5 



ergs 1), 



corresponding to FR I radio galaxies and most BL Lacs. These low jet power 'weak jet' sources, 
thought to exhibit radiatively inefficient accretion, are distinguished from the population of non- 
decelerating, low synchrotron-peaking (LSP) blazars and FR II radio galaxies ('strong' jets) which 
are thought to exhibit radiatively efficient accretion. The two-population interpretation explains the 
apparent contradiction of the existence of highly core-dominated, low-power blazars at both low and 
high synchrotron peak frequencies, and further implies that most intermediate synchrotron peak (ISP) 
sources arc not intermediate in intrinsic jet power between LSP and high synchrotron-peaking (HSP) 
sources, but are more misaligned versions of HSP sources with similar jet powers. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: jets 



1. INTRODUCTION 

The success of the current model for radio-loud Ac- 
tive Galactic Nuclei (AGN) as accreting super-massive 
black holes with bipolar, relativistic jets has rested on 
the ability of this physical paradigm to successfully ex- 
plain the chameleon nature of the spectral energy dis- 
tribution (SED) of this family of extreme objects. The 
jet model is inherently anisotropic and allows us to unify 
the highly luminous and variable blazars as relativisti- 
cally bea med counterparts of Fa naroff-Riley (FR) radio 
galaxies (jFanaroff fc Rilevlll974[) . 

The typical blazar spectrum exhibits a double peak 
(in vFu), attributed to synchrotron and inverse Compton 
(IC) radiation from the relativistic jet, which dominates 
in blazars due to the alignment of the jet axis near to our 
line-of-sight (f or a review of leptonic blazar models, see 
lB5ttcheill2007f ) . The synchrotron component peaks from 
sub-infrared (IR) energies up to hard X-rays, with the IC 
component peaking at much higher energies (typically 
gamma-rays). Spectroscopically, blazars are tradition- 
ally divided between optically 'featureless' BLLac ob- 
jects and (broad) emission line-emitting Flat Spectrum 
Radio Quasars (FSRQ). The long-standing oricntation- 
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based unification scheme for radio-loud AGN proposes 
BL Lacs and FSRQ as counterparts to FR I and FR 
II radio galaxies, respectively, based on similar spectra, 
morphologies, and rang e in extended radio luminosity 
(jUrrv fc Padovanil[l995h . 

However, violations to this scheme are well-known, 
and findings of p owerful, FR Il-lik e BL Lacs, and 
low-p ower FSRQ (iLandt et all 120061: iLandt fc Bignalll 
120081: iKharb et all I2O1O0 appear to break the sim- 
plest dichotomy (going back to the original FR classi- 
fication) between the lobe-dominated, high-power FR 
II galaxies and th e low-power, edge-darkened FR Is. 
IKharb et al.l ()2Q10[ ) also found BL Lacs displaying SEDs 
and hot-spots typical of quasars while others have ob- 
served BL Lacs exhibiting broad lines in low contin- 
uum states (iVermeulen et al.lll995t iCorbett et"aDl2000l : 
iRaiteriet al.ll2007f l. 

Further, the importance of optical spectral type in 
connecting radio galaxies with their aligned counter- 
parts is difficult to assess due to the strong effects 
of relativistic beaming on equivalent width measure- 
ments on which the blazar c lassifi cat ions are based 
(I Georganopoulos fc Marsched 119981 : iGhisellini et al.l 
120101 ). At the same time, the designation of Low and 
High Excitation Radio Galaxies (LERG and HERG, cor- 



Meyer et al. 



responding to prominence of e mission lines as measu red 
with an excitation index, see iHine fc Longaiii I1979D is 
increasingly found to be mixed among FR I and FR II 
morphological types; all of this creates serious problems 
for a unification scheme which connects spectral type 
with the morphology. 

1.1. The Blazar Sequence 

The synchrotron peak frequency (I'pcak) may take on 
a wide range of values (from 10^^ to 10^^ Hz at the ex- 
tremes), and is one of the principle ways to classify indi- 
vidual blazars. While FSRQ sources are typically found 
to have low i/peak (< 10^^-^ Hz), BL Lacs are found to 
span the entire range, with low, intermediate, and high- 
peaking BL Lacs (LBL, IBL, and HBL) typically dom- 
inating surveys selecting in the radio, optical, or X-ray, 
respectively. We adopt the generic terms for low, inter- 
mediate, and high synchrotron-peaking (LSP, ISP, HSP) 
blazars independen tly of the spectroscopic type, as in 
lAbdo et all (l2010d). 

iFossati et al.l ( 19981 ) found an anti-correlation between 
the frequency of the synchrotron peak and the bolomet- 
ric and synchrotron peak luininosities and suggested a 
link between the power of the source and fpeak- This 
sequence appeared to match a spectroscopic one, from 
predominantly powerful FSRQ sources at the low-peak, 
high-luminosity end th rough LBLs an d finally HBLs, at 
the low-luminosity end. IGhisellini et a l. (1998) suggested 
that more efficient cooling of particles in the jets of high 
luminosity blazars is responsible for the lower peak fre- 
quencies. However, it is not clear how the continuous 
blazar sequence fits in with the bi-modality of morphol- 
ogy and spectral type in radio galaxies. 

1.2. Beaming and Quantitative Unification 

Orientation-based unification predicts that misaligned 
blazars will drop in luminosity and synchrotron peak 
frequency according to the decrease in Doppler boost- 
ing with increasing angle. In the simplest case (a single 
Lorentz factor and a convex spectrum), the ratio of de- 
crease in log frequency to log luminosity will be approx- 
imately 1:4, roughly orthogonal to the sequence. (This 
holds for moving features in the jet; fo r a standing shock 
model the relation is 1:3, see e.g., iLind fc BlandfordI 
I1985D . The connection between blazars and the nuclear 
emission in radio galaxies can naturally inform our jet 
models. The synchrotron nature of the nuclear radio 
galaxy emission has been confirmed by comparison of 
the broad-band spectral indices and SEP characteris- 
tics of FR I nuclei with BL Lacs (IChiaberge et arill999l: 
Capctti e t "ani2000l: IChiaberge et al.ll2000HCapetti et all 
2002; Tru ssoni et al.ll2003[ ). where the effects of relativis- 
tic beaming (rather than obscuration) accounted for the 
10—10^ difference in radio and optical luminosities be- 
t ween the BL L a cs and radi o galaxies. 

iCapetti et all (|2002f ) and iTrussoni et all (|2003[) agree 
upon jet bulk Lorentz factors of F '^ 4—6 in order to 
account for the di fference in radio galaxy and blazar 
luminosities, while iHardcastle et al.l (120031 ) found F ^ 
3 sufficient to match the luminosity functions of the 
B2 and EMSS and 1 Jy samples. These values agree 
with VLBI studies of the pc-scale jets of HSP BL Lacs 
which show that they are barely super -luminal, with av- 
erage bulk Lorentz factors F ~ 2—3 ([Finer et al.l 1200 Jl : 
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2008ll2010D : however, se veral LSP BL Lacs s how high ve- 
lociti es on VLBI scal es (iJorstad et al.l 200lHCohen et ahl 
I2007D . Interestinglv. IHardcastle et all (|2003f ) also found 
that they could not reproduce the redshift distribution 
of 1 Jy BL Lac sample unless some of the high-redshift 
BL Lacs corresponded to FR lis. 

In contrast, high Doppler factors (~ 50 or greater) 
must be invoked to explain the brightness and rapid 
variability of TeV blazars, as the jet must be highly 
relativistic i n order to avoid absorpt i on by co-spatial 
IR photons (IDondi fc Ghiscllini 1995; Aharonian et al.l 
[20071: iGhiselhni et al...2009: .Giannios et al...2009. ). The 
discrepancy in apparent jet speeds and radio galaxy nu- 
clear emission can be resolved at least for FR I/BL 
Lac sources by invoking the presence of structured 
jets with multiple emitting regions of varying speeds. 
This inclu des the spine - sheat h axisymmetric m o del ex - 
plor ed bv iSwain et all ([19981) . IChiaberge eFall ([2000t) . 
a nd IGhisellini et al.l (I2005D . and the decelerating model 
of lGeorganopoulos fc KazanasI (|2003i ). all of which imply 
a high-speed spine or base which dominates the fiux when 
the source is aligned, and a slower component which dom- 
inates if the source is misaligned. Other blazar stud- 
ies also support the decelerating model for low power 
blazars, as recent VLBI observations of the MOJAVE 
samplc|3 reveal that the these sources also have lower ap- 
parent jet speeds, versus the high-power FSRQ which 
maintain coUima ted jets and consi stent jet speeds over 
longer distances ([Lister et al.ll2009( ). 

1.3. The Blazar Envelope 

The simplest scenario behind the appearance of a 
blazar sequence is that the physical properties of extra- 
galactic jets (as measured through the synchrotron peak 
frequency and luminosity) are a function of a single pa- 
rameter, the jet kinetic power (Lkin)- If there is a one-to- 
one correspondence between i^poak and Lkin, the original 
blazar sequence could be seen as the extreme 'edge' of 
an envelope, consisting of the most aligned (and there- 
fore brightest) objects. More misaligned sources should 
appear below the sequence, causing it to widen. If the 
jets are monoparametric in terms of jet power, these 
misaligned sources should fall in a line away from their 
aligned counterparts on the sequence, according to the 
drop in Doppler factor, forming 'tracks' of similar kinetic 
power in the envelope of blazars populating the I'poak - 
Lpoak plane. 

The original blazar sequence was based on 126 objects, 
a few dozen of which were detected in gamma-rays. Over 
a decade later there are thousands more blazars and 
blazar candidates in the literature and a much greater 
density of multi- wavelength data in general. Many of the 
newer samples have been designed to combat the early se- 
lection effects of the radio and X-r ay surveys, such as the 
low r adio power CLASS sample ([Caccianiga fc Marchal 
120041 ) ■ and it is now clear that as our surveys go deeper, 
many more low-luminosity blazars are emerging. These 
newer samples challenge both the sequence and its sim- 
plest extension that considers the sources in the enve- 

^ MOJAVE (Monitoring Of Jets in Active galactic nuclei with 
VLBA Experiments): Over 250 bright blazars being continuously 
monitored with VLBI. (www.physics.purdue.edu/mojave) 
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lope simply to be de-beamed bl azars. Contrary t o what 
is anticipated by the sequence, iLandt fc Bignalll (|2008[ ) 
showed that HSP and LSP BL La cs objects cover similar 
ranges in extended radio power. iCaccianiga fc Marchal 
(12004*) found several low luminosity, low t'pcak sources 
with relatively high core dominance compared to brighter 
objects at similar low i^peak, suggesting that the former 
are not simply related to the latter through misalign- 
ment, as one might naively expect. 

In this paper we examine the available evidence for 
a 'blazar envelope' beneath the upper-limit of the most 
aligned blazars, due primarily to mis-alignment. We test 
the possibility that jets are monoparametric engines by 
examining how the source power is related to other prop- 
erties of the blazar such as location in the fpeak — Lpoak 
plane, and re-evaluate the blazar sequence in light of 
the great increases in multi- wavelength, multi-epoch data 
available for blazar and radio galaxy jets. 

The paper is organized as follows: in ^ we explain 
the sample selection and describe our statistical and 
phenomenological models which allowed us to estimate 
the important envelope parameters for blazars and radio 
galaxies. In |J3] we present the results of the first round 
of data collection and analysis and the 'blazar envelope'. 
In g] we discuss the implications of the blazar envelope 
and the potential for constraining jet models. In iJ51 we 
summarize our findings. Throughout the paper we adopt 
a cosmology based on the concordance model with Qm 
= 0.27, f^A = 0.73, and Hq = 71 km s^^ Mpc^^. Energy 
spectral indices arc defined such that f^, oc i^^", and all 
luminosities (Li,) are references to powers (i^L^ implied). 

2. METHODS 

In order to investigate the nature of the blazar enve- 
lope, we require a sample of blazars with a well-sampled 
SED such that a reliable estimate of the synchrotron peak 
location can be made. The scarcity of such blazars may 
recommend utilizing as many sources as are known with- 
out regard for flux-limited samples, but this approach 
risks losing the ability to understand the limits of our ob- 
servations. We compromise by combining multiple well- 
defined (flux-limited) samples of flat-spectrum sources. 
The overall sample definition is thus complex, but does 
allow for the examination of biases by population mod- 
eling, as discussed in 



2.1. The Blazar Sample 

The flat spectrum candidate sample was made from the 
surveys listed in Table [T] A multi-frequency database of 
all sources in this master list was compiled by a combina- 
tion of all the data in the original publications, with ad- 
ditional data from the literature and from the SIMBAD 
and NED online databases, for a total of 3773 sources 
(taking into account duplicates). For 69 sources, previ- 
ously unpublished Swift data were available which were 
utilized to produce fluxes and X-ray indices. Spectro- 
scopic classification was taken from NED, SIMBAD, or 
the most recent publication, and non-blazar sources were 
removed, leaving 3014 blazars. A further reduction was 
made by removing sources with sparse multi-frequency 
SED sampling. The minimum required sampling was 
photometry at 3, 2, and 1 or more unique frequencies for 
the (rest frame) frequency ranges < 10^^'^, 10^^'^ — 10^^, 



and > 10^^ Hz, respectively, yielding 1327 sources suit- 
able for further analysis. We discuss the properties of 
the discarded sources in i j4.2.2l 

For sources with unknown rcdshift (~ 16% of the sam- 
ple) , we use the relationsh ip between optica l host galaxy 
magnitude and rcdshift of iSbarufatti et aD (j2005| ) to es- 
timate a value for these sources, taking the lower limit 
when exact host magnitudes were unavailable. All lu- 
minosities and frequencies were corrected to rest-frame 
values before spectral fitting; the global fits (rather than 
an assumed index) were then used to find k-corrected 
rest- frame quantities of interest. 

2.2. Jet Power 

Measuring the jet kinetic jet power of our blazars and 
radio galaxies allows us to see if the areas occupied by jets 
of similar power in the z/pcak ~ Lpoak plane can be under- 
stood simply as an orientation effect (as it would be an- 
ticipated in the original blazar sequence) , or if additional 
considerations are needed. We note that for our purposes 
we need only classify sources into a few broad bins of jet 
power (as we show below, a width of one decade gives 
us a sufficient number of well populated jet power bins), 
and an extremely accurate estimation of the jet kinetic 
p ower is not needed. 

iRawlings fc Saunder j ()199H) estimated jet power by di- 
viding the lobe energy content of radio galaxies by their 
age. The energy content was estimated through mini- 
mum energy arguments {e.g. iMilevI [l980[ ) derived from 
measurements of the low frequency radio emission where 
the lobe emission is assumed to dominate over the jet, 
while the source age was estimated from spectral breaks 
in the lobe spectra th rough synchrotron aging arguments. 
iWillott et al.l (|1999[ ) revisited the issue of obtaining the 
jet power from the low frequency extended luminosity. 
They studied in detail the parameters that come into this 
relation, including the unknown energy content of the 
non-radiating partic les and the time evolu tion of the lobe 
luminosity (see also iBlundell et al]|1999[ ) and presented 
a scaling relation which included a careful parametriza- 
tion of the various uncertainties (see their Section 4.1, 
and Equation 11). 

The actual normalization for this scaling between the 
jet kinetic power and the low frequency extended lumi- 
nosity (Lcxt) came from a very different route. The dark 
cavities seen in (e.g., Chandra) X-ray images are created 
by AGN jet activity and can be used to estimate the pAV 
work needed to inflate them. Importantly, this energy 
estimate does not depend on the minimum energy as- 
sumption and does not require knowledge of the amount 
of non-radiating matter in the lobes. The jet power Lkin 
is estimated by dividing the pAV work needed to ex- 
pand the cavity by t he cavity age found by estimating the 
lobe buoyancy time (jBirzan et al.ll2008t iCavagnolo et al.l 
[2OT0h . 

While this method is limited to a small number of 
sources at present, the Lkin — Lcxt relation covers al- 
most six orders of magnitude in jet power, including 
both FR I and FR II sources. The slope of the scaling 
ICavagnolo et al.l (l2010f) flnd is co r npatib le within errors 
with that found bv IWillott et al.l (|1999[ ) and requires a 
lobe energy content larger by at least a factor of 100 com- 
pared to that inferred by minimum energy arguments. 
The cavity work shows convincingly that Lkin and the ex- 
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TABLE 1 
Blazar Samples 



Sample 



Abbrev. No. of 

Objectst 



References 



1 Jansky blazar sample 
2-degree field (2dF) QSO survey 

2 Jansky survey of flat-spectrum sources 
The Candidate Gamma-Ray Blazar Survey 
Cosmic Lens All-sky Survey 

Deep X-ray Radio Blazar Survey 

EGRET blazar sample 

Einstein Slew survey sample of BL Lac objects 

Fermi blazar sample 

Hamburg-RASS bright X-ray AGN sample 

Parkes quartcr-Jansky flat-spectrum sample 

Radio-Emitting X-ray source survey 

Radio-Optical-X-ray catalog 

RASS - Green Bank BL Lac sample 

Einstein Medium-Sensitivity Survey of BL Lacs 

RASS - SDSS flat-spectrum sample 

Sedentary survey of high-peak BL Lacs 



1 Jy 


34 


2QZ 


56 


2 Jy 


213 


CGRaBS 


1625 


CLASS 


236 


DXRBS 


283 




555 




48 


IFGL 


689 


HRX 


104 




878 


REX 


55 


ROXA 


816 


RGB 


127 


EMSS 


44 




501 




150 



Stickel et aP lfl99il) 



Londish et al.l 112002!. [^OOTD 

SVall & Peacock ( 1983) 

Hcalcv et al. (2008) 

Marcha et al. ( 20011') : ICaccianiga eFaLl iWO^) : 

Caccianiga & Marchal p 004l ) 

Perlman et al. (199?); ILandt et all | |200!1) : 

Padovani ct al. (2007J 

Mukherjce ct al. ( 1993) 

Perlman et al. ( 1996) 

'Abdo ct al. (2010a) 

[Beckmami et a l. (20031) 

IJackson et all (120021 ): JHook et al.l | |2003) : IWall et al.l 

J2005I) 

(Jaccianiga ct al. (19991) 
furriziani et al. (200 3) 
Laurent-Muehlciscn ct al.l (119991) ^_^_ 
Morris ct al. (1991); Rec tor et al.l | |2000|) 
Plotkin et al.i (i2008i) 



Giommi et"all ( IT999I. [20051) : [Piranomonte et~all ( 120071) 



t Number of sources listed is total sources in sample (not necessarily number of blazars) . 



tended radio luminosity are related. We therefore chose 
the low-frequency extended luminosity (also commonly 
referred to as steep or lobe emission) at 300 MHz (L300) 
as an estimator of the jet power for this work. This value 
is then converted to Lkin from the linear fit (with l-cr er- 
ror), 



log Lkin = 0.64(±0.09) (log L 



300 



40) 



-t-43.54(±0.12) [ergss"^] 



(1) 



taken from lCavagnolo et al.l ()2010[ ). Although the scatter 
in the scaling is not insignificant (^ 0.7 decade in Lkin), it 
is adequate for our purpose of classifying sources in broad 
bins of jet power. This value (Lkin) is what we generally 
refer to throughout the paper as the jet (kinetic) power, 
not to be confused with Lpeak or other observables of the 
jet which are affected by beaming. 

2.2.1. Extended Radio Emission 

The extended radio emission from the lobes of radio 
galaxies differs from the beamed core emission in be- 
ing both non-variable and generally exhibiting a much 
steeper spectrum. Thus one way to isolate the extended 
emission is to look for a change in the spectral index 
at low frequencies. For our sample, a spectral decom- 
position (SD) method was applied to all sources which 
had data points at five or more frequencies below 10^^ 
Hz. Sources were tested for a two-component spectrum 
by fitting alternately with a double power-law and then 
a single power law (assuming normal error). Because 
the latter can be expressed as a particular case of the 
former, this allowed us to use the maximum likelihood 
(L) ratio test for nested models. In this test, the statis- 
tic — 2Alog(L) distribution can be approxim ated as a x^ 
with 712 - ni degrees of freedom (jWilkal 19381 ) . The double 
power-law model has 712 = 4 degrees of freedom (slope 
and normalization for each powerlaw) , while the null hy- 
pothesis model has ?ii = 2, for a total of two degrees of 



freedom. We rejected the two power-law fit for a signif- 
icance level a >0.05. The best steep-spectrum fit from 
the hypothesis test was then used to calculate L3oo- 

The normal approximation for the spread around the 
power law is a rough approximation, and for eight sources 
it appears the fit failed to find a significant two-powerlaw 
solution due to a great deal of variability in the radio 
core points. L300 was estimated for these sources by 
fitting hand-selected low-frequency radio points corre- 
sponding to a stable emission component persistent be- 
neath the core variability. This method was also used 
for a few sources which passed the test, where core vari- 
ability caused the steep power law fit to be somewhat 
more fiat than would be fit by eye (a correction of a few 
percent). We note that inclusion of these eight sources 
and the corrections do not materially affect the results 
presented. 

The SD method is the most reliable way to separate the 
boosted core emission from the isotropic lobe emission, 
as it is model independent. Another method is to mea- 
sure the extended emission from core-subtracted VLA 
images (typically at 1.4 GHz). In order to expand our 
sample and counter some selection bias (see longer dis- 
cussion in ij4.2p . extended flux measurements were taken 
from the literature where available (see Table [5] for ref- 
erences). One difficulty with these measurements is that 
an assumption must be made about the spectral index 
used to extrapolate from 1.4 GHz to the 300 MHz lumi- 
nosity. It is also clear from sources with both a visible 
steep spectrum and VLA data that these are not always 
in agreement. Even for the same object, VLA estimates 
by different authors can vary widely, as in the case of 3C 
446, with 1.4 G Hz extended flux estimates of 91.6, 271, 
and 4030 mJv (iKharb et all [Mot iCooper et all 120071: 
lAntonucci fc UlvestadI 119851 )7 none of which agrees with 
a clear steep component which suggests a flux ~ 800 mJy 
at that frequency. However, examination of the over 100 
sources with both VLA and SD estimates reveals that 
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Fig. 1. — SEDs for three sources in which the extended radio 
emission is clearly visible as a steep component at low frequen- 
cies. The FSRQ QSO B0923+392 (top) has a very clear steep 
component at low frequencies, dominating up to ~ 10^'^ Hz, and a 
prominent 'big blue bump' from the accretion disk, suggesting that 
it may be relatively misaligned compared with the Fermi-detected 
BL Lac IJy 2007-1-77 (middle), because they differ by ~ 2 decades 
in extended luminosity but only ~ 1 decade in beamed power. This 
latter source has a similar extended radio spectrum to the radio 
galaxy 3C 338 (bottom), however in the radio galaxy more of the 
extended emission is visible due to the nuclear emission being much 
less beamed. We also show estimates of the extended emission at 
1.4 GHz as measured from VLA maps (diamonds), which generally 
agree with the spectral fits. The dashed curves are the best-fit of 
the blazar model described in the text and in Appendix |A] 

typical differences are less than 0.5 in log luminosity, with 
the best agreement found using a low-frequency index of 
a = 1.2 for the extrapolation to 300 MHz. 

The full SED for three sources is shown in Figurc[TJ the 
bright quasar QSO B0923+392 (top), the LSP BL Lac 
IJy 2007+77 (middle), and the FR I radio galaxy 3C 
338 (bottom). In all cases, extended radio emission was 
detected in our two-component SD test - these points 
are shown in lighter colors of pink, green, and orange, 
respectively. (Estimates of the extended luminosity de- 
rived from VLA maps are also shown as diamonds). We 
note that while IJy 2007+77 and 3C 338 have similar ex- 
tended luminosities at 300 MHz, the much greater core- 
dominance of the blazar is evident, and results in much 
less of the extended emission being visible. Also, while 
the quasar is evidently more powerful than the other two 
sources by two orders of magnitude in L300, we see that 
it is likely more misaligned than IJy 2007+77, from the 
greater visibility of the steep component and the promi- 
nent accretion disk bump in the optical. This relatively 
bright source was also not a significant detection in the 
Fermi/LAT 1-year catalog, unlike IJy 2007+77. This 
could be explained if the gamma-ray emission is more 
beamed than that at lower frequencies. 

2.2.2. Connection Between Extended and Unheamed Core 

Emission 

For the SD method, it is important to note that failure 
to find a two power-law solution for a source can just as 
easily be due to a lack of low-frequency data coverage as 
to the source emission being highly jet dominated down 
to low frequencies. The approximate frequency at which 
the beamed emission becomes dominant varies from 10^ 
to lO^'' Hz in our sample. In Figure [2] we plot the 'cross- 
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Fig. 2. — The crossing frequency, i^cross is the approximate fre- 
quency at which the SED becomes dominated by the beamed jet 
emission rather than the steep extended emission. Radio galaxies 
are shown as empty squares, and blazars as filled circles. As log 
core dominance (Rce) increases due to beaming this crossing point 
moves to lower energies. Many sources will remain core-dominated 
to low frequencies, even to the limit of ground-based observations 
at ~ 10 MHz. The broken line shown is the linear correlation be- 
tween the plotted variables (r = 0.87). The solid red line marks the 
300 MHz frequency at which the extended emission is measured. 

ing' frequency i^cross of the fitted jet SED and steep com- 
ponent and the log core dominance Rqe (we define Rce 
as the log of the ratio of core to extended luminosity at 
1.4 GHz). The correlation is well-fit by a linear regression 
model (correlation coefficient r = 0.87). This correlation 
can arise from two coexisting causes: an intrinsic spread 
in unbeamed core power for a given extended luminosity 
and the effect of beaming. If the effect of the natural 
spread of core power dominates over the effects of beam- 
ing (as would be the case if beaming of the core emission 
is weak or absent), then we expect to see radio galaxies 
and blazars mixed along this correlation (this is because 
the range by which each particular source moves before 
being classified as a radio galaxy will be smaller than the 
initial intrinsic range that the aligned blazars populate). 
This is, however, not what we find, as shown in Figure 
[2I Instead, we observe a clear separation of radio galax- 
ies and blazars which suggests that the effect of beaming 
dominates over the intrinsic spread and that the range of 
unbeamed core emission for a given extended luminosity 
be narrow. 

2.2.3. Caveats 

Many of our sources which lack VLA imaging have ob- 
servations down to very low frequencies (30 — 200 MHz), 
and show no signs of a steep component. It is impossible 
to tell if these sources are highly core-dominated to the 
lowest observed frequencies or if the lobe emission has 
a flat spectrum which is indistinguishable from the core 
emission. The latter scenario may be a concern for the 
very low power sources, as it is possible that the break 
frequency (i.e., the point at which the lobe emission be- 
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Fig. 3. — Examples of typical blazar SEDs in the master sam- 
ple of flat-spectrum sources. Visual inspection of over 1300 SEDs 
was used to separate sources with unambiguous synchrotron peaks 
(top, 653 sources total) from those with under-sampled SEDs which 
might be consistent with several possible peak locations (bottom). 
The dark line is the best-fit model (described in Appendix |A]| for 
each case, while the dashed lines show other possible SED shapes 
consistent with the more sparsely-sampled SED. 

comes st eep) may increase wi th decreasing synchrotron 
power as lBrunetti et al.l (J2003D found for hot spots in ra- 
dio galaxies. For these sources it may be advisable to 
perform VLA map subtractions at multiple frequencies 
to confirm such a spectrum. 

Finally, we must point out the risk in the common 
practice of assuming that flux below ^ 300 MHz, taken 
without examination of the spectral characteristics, is 
actually from the extended isotropic component, rather 
than the core. This assumption would be clearly inap- 
propriate for many of our sources, as shown by Figure [21 
and could lead both to an overestimate of the jet power 
and an underestimate of the radio core dominance often 
used to infer the degree of mis-orientation of the source. 

2.3. Fitting the Synchrotron Spectrum 

For the 1327 blazars making the initial cut based on 
SED sampling, we utilized a phenomenological paramet- 
ric model in order to estimate the synchrotron peak fre- 
quency and luminosity. The basic shape is parabolic, 
with a power-law extension to low frequencies and a ris- 
ing power- law to sim ulate the begiiining of the IC spec- 
trum in the X-rays (jFossati et al.l I1997D . These com- 
ponents were free to vary within a reasonable range in 
shape or slope, respectively, and in relative normaliza- 
tion, for a total of 6 parameters. In sources where a 
big blue bump was suggested on visual inspection, this 
feature was added to the SE E) fit by adding a template 
accretion-disk spectrum (see iFrank et al.l 120021 ). adding 
one additional parameter for the bump peak luminosity 
(see Appendix |A] for a summary of the model) . 

To fit our model, we used a simulated annea ling rou- 
tine for the parameter search (jGofFe et al.l 1 1994 ). and an 
algorithm based on a combination of minimum square 
distance (in log L^) between model and data, with the 
added constraint of reasonable agreement with the X-ray 



spectral index where available, enforced by increasing the 
score of the fit by a function of the angle between the 
model and observed spectral slope. 

We note that while simultaneous multi- wavelength ob- 
servations give the best 'snapshot' view of an entire SED, 
such data is usually still biased by being limited to a 
single epoch, as the variability of these objects is well- 
known. The compilation of all available photometry, gen- 
erally collected over a long time-frame, sampling many 
states of the SED, allows us to fit the model to the av- 
erage SED of the source which may be more appropriate 
in a population context, in which extreme source states 
might confuse results. 

2.4. The TEX Sample of Blazars 

After fitting with the SED model, each source was sub- 
sequently examined by eye, in order to identify ambigu- 
ous cases where the formal best-fit model did not clearly 
represent the only reasonable SED fit. In many cases 
these removed sources had few data points in the critical 
region from IR to UV, lacked X-ray spectral indices, or 
had significant host galaxy contamination in the optical. 
A comparison of an acceptable and a rejected source with 
their best SED fits is shown in Figure [3l A total of 653 
sources remained after this examination, in which we feel 
confident in the determination of the synchrotron peak 
to within one order of magnitude. (The rejected sources 
are discussed along with the other under-sampled blazars 
in §4.2.2p . In total, 216 sources have robust estimates 
of t'pcak, Lpcak, and L300; this set comprises the TEX 
(trusted extended) sample. Of these, 129 have extended 
radio fiux from spectral decomposition, 87 were added 
using only VLA map data, and 90 have both a spec- 
tral and VLA estimate. For the remaining 437 blazars 
comprising the UEX (unknown extended) sample, upper 
limits on L300 (and subsequently Lkin) were found by ex- 
trapolating from the luminosity at the lowest available 
frequency using a = 1.2. 

In Table [2] we list the sources which have a con- 
firmed synchrotron peak and extended power measure- 
ment (TEX sample). In columns 1 — 5 wc give the 
source name, position, optical type, and redshift. Es- 
timated redshifts are noted by an asterisk in column 5. 
In columns 6 and 7 we give the best-fit synchrotron peak 
frequency and luminosity, respectively. In columns 8 and 
9 we give L300 and the converted value of Lkin, respec- 
tively. In column 10 we note the method for the extended 
radio flux estimate: VLA, spectral decomposition (SD), 
or both. We note those sources (19 total) in which the 
spectral result is used instead of a VLA measurement by 
a dagger in that column. In column 11 wc give the ref- 
erences for VLA data. All values are in the rest-frame of 
the source. 

Because of the somewhat subjective requirement of an 
unambiguous synchrotron peak, the TEX sample does 
not contain all the sources in original list, and cannot be 
considered a complete sample. We focus most of our dis- 
cussion on qualitative analysis and positive flndings for 
our well-characterized sample. Relaxing the criteria for 
inclusion in the sample risks inaccurate measurements of 
the three principle quantities (i^peak, Lpcak, and Lkin) and 
could lead to incorrect conclusions. See t j4.2l for a discus- 
sion of the selection effects of our master sample, the 
nature of the rejected sources, and the possible impact 
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TABLE 2 

TEX Combined Blazar, Source List 



Name 


R.A. 


Dcd. 


Opt.t 


z 


log I^pcak 


log Lpeak 


log L300 


log Lkin 


Method 


Ref. 




(J2000) 


(J2000) 






(Hz) 


(ergs s-i) 


(ergs s-i) 


(ergs s-i) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


QSO B0003-06 


00 06 13.89 


-06 23 35.33 


B 


0.347 


13.16 


46.22 


41.68 


44.21 


both 


f,h 


QSO B0016+73 


00 19 45.73 


+73 27 30.07 


Q 


1.781 


13.26 


47.65 


42.61 


44.74 


both 


f 


QSO B0048-09 


00 50 41.31 


-09 29 05.21 


B 


0.200 


13.43 


44.97 


41.20 


43.93 


both 


b,d,f 


4C 01.02 


01 08 38.77 


+01 35 00.27 


Q 


2.099 


13.00 


47.71 


44.57 


45.85 


both 


c,f,h 


QSO B0109+224 


01 12 05.82 


+22 44 38.78 


B 


0.265 


13.87 


45.57 


40.33 


43.44 


both 


b,f,h 


QSO B0118-272 


01 20 31.66 


-27 01 24.65 


B 


0.557 


14.04 


46.13 


42.38 


44.61 


VLA 


d 


QSO B0133+47 


01 36 58.58 


+47 51 29.30 


Q 


0.859 


13.34 


46.63 


41.93 


44.35 


VLA 


f,h 


QSO B0138-097 


01 41 25.83 


-09 28 43.67 


B 


0.733 


13.49 


46.22 


42.41 


44.62 


VLA 


d 


4C 15.05 


02 04 50.41 


+ 15 14 11.04 


Q 


0.405 


12.39 


45.15 


41.39 


44.04 


VLA 


h 


QSO B0212+73 


02 17 30.82 


+73 49 32.55 


B 


2.367 


13.38 


47.20 


42.30 


44.56 


VLA 


f,h 


QSO B0215+015 


02 17 48.95 


+01 44 49.69 


Q 


1.715 


13.72 


47.71 


43.52 


45.26 


VLA 


b,f,h 


QSO B0218+357 


02 21 05.48 


+35 56 13.70 


Q 


0.960 


12.86 


45.91 


43.09 


45.01 


SD 




3C 66A 


02 22 39.61 


+43 02 07.80 


B 


0.444 


13.98 


46.34 


42.85 


44.87 


both 


b,e 



Note. — Table [2] is published in its entirety in the eleetronie edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 

Rcdshift estimated from optical magnitude 
t Spectral decomposition estimate of extended flux used instead of VLA estimate 
t Optical Spectral Type . B = BL Lac, Q = FSRQ. 
fa lAntonucci fc UlvestadI l lT985f) 
^ IMurphvet al.l II1993I) 
' ICassaro eTa P fidW) 
"^ ILandt fc B ignall (20Ql) 
^ IKharbet"aI.i(i20ia) 
E ICondon eTal.! I I2002I) 

'^ ICooper et"an ||2007D 

' ICaccianiga fc Marctial I I2004I) 



on our conclusions. 

2.5. Radio Galaxies 

As the jets isolated in radio galaxies are understood 
to be the less-beamed counterparts to blazars, one 
way to understand the relationship of the synchrotron 
peak and beaming in the envelope is to see where the 
synchrotron peak falls for these jets. In the context of 
the presumed blazar sequence, these 'hidden blazars' 
should indicate whether the sequence is preserved in 
misaligned jets. We take advantage of recent work to 
isolate the radio, IR, optical, and X-ray core emission 
from the jet in radio galaxies using high-resolution 
mapping to place as many rad io galaxies as possible 
on the i^peak - L pcak plane (IChiaberge et al.l 119991: 
Capetti et al.ll200(3l: IChiaberge et al.ll200a iCapetti et alj 



2002t iHardcastleet all 120031 'Trussoni ct al.' '2005 
Balmaverde et all l2006t [^lassaro et al. 201(| 
Buttiglione et al.l 120111 ). We examine the 45 radio 



Fourteen radio galaxies are nearby, well-known sources, 
with core measurements at multiple frequencies which al- 
low fitting of the synchrotron SED with the same param- 
eterization as for the blazars. This s et includes seven o f 
the Fermi-detected misaligned AGN (jAbdo et al.ll2010b[) . 

For the remaining 31 with sparser SED sampling, the 
synchrotron z/poak was estimated using a non-parametric 
likelihood estimator. The joint distribution of Rce, 
t'pcak, and the SED color^ Uro, ctox was calculated for 
all the sources in our TEX, UEX, and fitted radio galaxy 
samples. From this, a conditional, one-dimensional den- 
sity on t'peak can be found by supplying the other three 
observables. The value of the peak was taken as the 
maximum of this distribution, and we also report the 
20% and 80% quartiles as the sampling error. 

It is straightforward to estimate the peak luminosity 
from the radio luminosity at 5 GHz, with which it has a 
linear correlation (correlation coefficient r = 0.85). The 
linear fit (with 1-cr error) is 



galaxies with extended radio flux measurements, and at 
least one radio, optical, and X-ray measurement of the 
nuclear emission available in the literature. The L300 for 
all radio galaxies was found from a parabolic fii0 to the 
steep, low-frequency radio emission and subsequently 
converted to Lkin as in Equation [ij In all cases where 
available, VLA map-derived estimates at 1.4 GHz were 
in agreement with the spectral fits. 

^ For fitting the steep component in radio galaxies, a log- 
parabolic model is more appropriate than a power-law due to the 
much greater range of the spectrum visible, which usually has some 
curvature in uF,y. 



log Lpeak = 0.61 (±0.01) (log LsGHz " 43) 

+ 45.68 ±(0.02) [ergss" 



(2) 



The 90% confidence interval, assuming a normal er- 
ror distribution, corresponds to an uncertainty of A (log 
Lpcak) — 0.8. 

The most likely peak frequency found by our statisti- 
cal model for all radio galaxies is at or below 10^"* Hz, 
with the exception of four FR II sources with moder- 

^ Spectral indices are defined at 1.4 GHz, 5000 A, and 1 keV. 
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TABLE 3 
















Radio Galaxies 








Name 


R.A. 


Decl. 


Morph.t 


z 


log fpcak 


log Lpcak 


log L300 


Lkin 




(J2000) 


(J2000) 






(Hz) 


(ergs s-i) 


(ergs s-i) 


(ergs s-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C 380+ 


18 29 31.77 


+48 44 46.17 


2 


0.69 


12.95 


45.68 


44.06 


45.56 


3C 17 


00 38 20.52 


-02 07 40.8 


2 


0.220 


13.05«-j0 


45.35 


42.82 


44.86 


3C 18 


00 40 50.35 


+10 03 23.14 


2 


0.188 


12.83+°g 


44.79 


42.55 


44.70 


3C 29 


00 57 34.92 


-01 23 27.89 


1 


0.045 


12.94+0-39 


44.01 


41.34 


44.02 


NGC 315 


00 57 48.88 


+30 21 08.81 


1 


0.016 


13.96l^:'65 


43.94 


39.78 


43.13 


3C 31 


01 07 24.93 


+32 24 45.20 


1 


0.017 


13.09«i? 


43.53 


40.79 


43.70 


3C 33.1 


01 09 44.26 


+73 11 57.17 


2 


0.181 


12.99+111 


44.30 


42.42 


44.63 


3C 66Bt 


02 23 11.41 


+42 59 31.38 


1 


0.02 


13.36 


42.09 


40.73 


43.67 


3C 78t 


03 08 26.22 


+04 06 39.30 


1 


0.03 


12.99 


44.02 


40.94 


43.79 


30 83. IB 


03 18 15.73 


+41 51 27.38 


1 


0.025 


13.09l°;g 


43.40 


40.93 


43.78 


30 84t 


03 19 48.16 


+41 30 42.10 


1 


0.02 


13.46 


43.09 


40.16 


43.34 


30 lilt 


04 18 21.27 


+38 01 35.80 


2 


0.05 


13.28 


43.29 


41.41 


44.05 


30 133 


05 02 58.50 


+25 16 24.0 


2 


0.278 


13.04l0;i 


45.08 


42.98 


44.95 


30 135 


05 14 08.34 


+00 56 31.63 


2 


0.127 


13.05l°:l 


43.79 


42.04 


44.41 


30 165 


06 43 06.70 


+23 19 00.3 


2 


0.296 


13.04l?:^3 


44.33 


42.80 


44.84 


30 171 


06 55 14.87 


+54 08 59.4 


2 


0.238 


13.04l°;i 


43.89 


42.73 


44.81 


30 189t 


07 58 28.10 


+37 47 11.80 


1 


0.04 


13.44 


42.34 


40.88 


43.75 


30 264 


11 45 05.00 


+19 36 22.74 


1 


0.022 


12.60+°l2 


43.77 


40.72 


43.66 


30 270 


12 19 23.21 


+05 49 29.69 


1 


0.007 


13.W+Jil 


43.33 


40.09 


43.30 


NGC 4278 


12 20 06.82 


+29 16 50.71 


1 


0.002 


12.96lO;f7 


42.55 


37.29 


41.70 


M 84t 


12 25 03.74 


+12 53 13.13 


1 


0.00 


12.98 


41.19 


39.04 


42.70 


M 87t 


12 30 49.42 


+12 23 28.04 


1 


0.00 


13.01 


41.58 


40.95 


43.79 


Cen At 


13 25 27.61 


-43 01 08.80 


1 


0.00 


12.72 


41.47 


40.21 


43.37 


30 287.1 


13 32 53.27 


+02 00 45.79 


2 


0.216 


13.05+_rol 


45.21 


42.38 


44.61 


30 296t 


14 16 52.92 


+10 48 26.43 


1 


0.02 


12.99 


41.24 


40.57 


43.57 


30 300 


14 23 01.00 


+19 35 17.0 


2 


0.272 


-,r,7.,+0.40 


44.32 


42.85 


44.87 


30 317t 


15 16 44.48 


+07 01 18.00 


1 


0.03 


12.83 


42.95 


41.34 


44.01 


30 321 


15 31 43.47 


+24 04 18.93 


2 


0.096 


io 97+0.60 
-'-•^■^'_0.48 


44.19 


41.86 


44.31 


B2 1553+24 


15 56 03.87 


+24 26 52.70 


1 


0.043 


13.1ll??l 


43.79 


40.48 


43.52 


30 338t 


16 28 38.27 


+39 33 04.97 


1 


0.03 


13.12 


41.44 


41.18 


43.92 


NGC 625lt 


16 32 31.96 


+82 32 16.39 


1 


0.02 


12.79 


43.25 


40.68 


43.64 


30 346t 


16 43 48.62 


+17 15 48.97 


1 


0.16 


13.46 


43.49 


42.25 


44.53 


Her A 


16 51 08.14 


+04 59 33.32 


2 


0.155 


13.04l0-i 


44.08 


43.60 


45.30 


30 349 


16 59 32.00 


+47 02 16.0 


2 


0.205 


12.94i0;25 


44.39 


42.49 


44.67 


30 388 


18 44 02.41 


+45 33 29.81 


2 


0.091 


i3mt°:tt 


44.32 


42.00 


44.39 


30 442A 


22 14 45.00 


+13 50 47.66 


1 


0.026 


i2.98i°;i 


42.74 


40.77 


43.69 


30 449 


22 31 20.62 


+39 21 29.81 


1 


0.017 


13.1810J1 


43.27 


40.17 


43.35 


B2 2236+35 


22 38 29.46 


+35 19 40.01 


1 


0.027 


13.04l°;23 


43.03 


39.62 


43.03 


30 452 


22 45 48.78 


+39 41 15.36 


2 


0.081 


ism+fll 


44.48 


43.34 


45.15 


30 460 


23 21 28.50 


+23 46 48.0 


2 


0.268 


13.04i;igl 


44.52 


42.52 


44.69 


30 465 


23 38 29.39 


+27 01 53.52 


1 


0.031 


13.07t°oil 


44.07 


41.12 


43.89 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcmin- 
utcs, and arcscconds. 

t Synchrotron peak estimated from blazar model fit to core fluxes 
t Morphological Type. 1 = FR I, 2 = FR II. 



ate jet powers (Lkin ~ lO''^-^ ergs s~^) which display 
extreme X-ray luminosities. These sources include the 
broad-lined radio galaxies (BLRG) 3C 390.3 (which has 
been previously reported to have an LBL-like SED by 
I Ghosh fc Soundararaiaperumall |1995() and 3C 332, as 
well as two lesser-studied narrow line radio galaxies, 3C 
197.1 and 3C 184.1. The two BLRG have been de- 
tected with Swift and BeppoSAX in the hard X-ray 
(jGrandi et all 120061: iCusumano et al.l I2010D . Ah four 
sources were excluded from the final sample, as the X-ray 



is likely either thermal or from a strong IC component, 
and without knowledge of the synchrotron X-ray emis- 
sion, any peak estimate would be practically arbitrary. 

The final list of 41 radio galaxies is given in Table |31 
Column descriptions are as in Table [2] 

3. RESULTS 

3.1. The Blazar Envelope 

The TEX sample of 216 blazars, and the 41 radio 
galaxies are shown in the plane of Lpeak versus ;/peak in 
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Fig. 4. — The blazar sequence, which originally showed an anti- 
correlation between synchrotron Lpcak and i^pcak has been ex- 
panded into an 'envelope' with the addition of new observations 
and radio galaxies. BL Lacs are shown as filled circles, FSRQ as 
filled triangles, and radio galaxies as squares (FR I) and inverted 
triangles (FR II). Color indicates the jet kinetic power (L]ji„ in ergs 
s~^), as estimated from extended radio fiux measurements. Track 
(A) shows the path of a synchrotron peak for a single-component 
jet, and (B) for a decelerating jet of the type hypothesized to ex- 
ist in FR I sources. The fully-aligned limit for each (0 degrees) is 
shown as marked, with the arrow direction indicating the move- 
ment of the model source as it is misaligned. 

Figure |31 Blazars are represented as circles (BL Lacs) 
and triangles (FSRQ). The power of the source (Lkin) is 
shown by color in bins of one decade. Examining the 
location of blazars, there is an obvious lack of sources 
with both high jet powers and synchrotron peak fre- 
quencies above 10^^ Hz, though a few of the ISP/HSP 
sources have quasar-like spectra. The population at low 
peak frequencies, alternately, is highly mixed in terms of 
kinetic power (with many low-power sources coincident 
with high-power sources), spectral type (with both BL 
Lac and quasar spectra present) and the spread of peak 
luminosities, which ranges from 10'*'* to 10^^'^ ergs s~*. 
The 41 radio galaxies are shown in Figure |4] as squares 
(FR I) an d inverted triangles (FR II). Those detected 
by Fermi (jAbdo et al.l l2010bD are circled in red. For 
our radio galaxy sample, we note that the jet power 
does not apparently exert much influence on their syn- 
chrotron peak frequency, and it appears that all radio 
galaxies in our sample exhibit jets with peaks below the 
optical-UV range in frequency. Further, no evidence is 
found for any conservation of a 'blazar sequence' within 
the radio galaxy jets shown here. If, as has been sug- 
gested, different emission regions in the less powerful jets 



dominate at different orientation angles (iChiaberg e et al. 



2OO0 ; iGeorganopoulos fc Kazanail2003t iGhisellin i et al. 



2005f) it may simply be that the 'sequence' behavior (in 
which jet power regulates the synchrotron peak) pertains 
to the fast parts of the jet which are no longer visible in 
highly misaligned radio galaxies. 

It becomes interesting to ask the question of what 
type of orientation tracks are expected on the Lpcak- 
fpeak plane if we assume that powerful sources do not 



decelerate, while weak sources require velocity gradients. 
The former assumption is supported by the fact that 
the Lorentz factors required to model the blazar SED 
of these sources (F ~ 10 — 20) are similar with those 
dedu ced from following V LBI superluminal components 
(e.g. iJorstad et al.l 120011 ) . while the latter (as we dis- 
cussed in ijl.2p is a result of our need to (i) explain the 
discrepancy between the high Lorentz factors required 
by the fast variability, pair production opacity, and SED 
modeling for the sub-pc scale jet and the lack of signifi- 
cantly superluminal motions at the VLBI scales and (ii) 
satisfy the unification between low power blazars and FR 
I radio galaxies. 

Two tracks are shown in Figure ID These illustrative 
tracks do not intend to simulate any particular source, 
and are shown to simply demonstrate how sources that 
start with the jet pointing toward us will appear as the 
jet orientation angle increases. The first, labeled (A), 
follows the misalignment of a powerful source assuming 
a single jet Lorentz factor F = 15. This track follows 
the expected 1:4 ratio of decreasing (log) frequency to 
(log) luminosity of the p eak as the source is decreasingly 
beamed (see, e.g.. iLind fc Blandford 198f|). The second 
track, labeled (B), is for a weak source and shows the 
change in peak location for one possible model of a decel- 
eratin g jet, as described in IGeorganopoulos fc Kazanai 
()2003[) . In this particular example, the Lorentz factor 
of the flow decreases continuously from T„iax = 15 to 
Fmin = 3 in a distance of 5 x lO*'' cm, and the emitted 
spectrum is calculated along the jet and integrated spa- 
tially. The track exhibits the qualitative characteristics 
{e.g., more horizontal movement) not only consequent 
of the decelerating model, but of any model with veloc- 
ity gradients in which the most energetic electrons are 
injected in the fast part of the flow, s uch as the spine- 
sheath model of IGhisellini et all (|2005| ). Such models, as 
can be seen in Figure JH exhibit a shallow decrease of 
Lpcak and a sharp decrease of i^peak because the opening 
angle of the beaming cone is smaller for higher frequency 
emission; as the angle to the line of sight increases, pro- 
gressively lower frequency emission is beamed away from 
the observer. Only for angles > 1/F„„;„, do we expect 
the beaming track to start resembling that of a constant 
Lorentz factor flow and Lpoak to start dropping faster 
with decreasing t'peak- This model predicts that the low- 
power HSP blazars will follow a more horizontal path 
from the most beamed blazars to the misaligned radio 
galaxies. Interestingly, though each model originates in 
very different parts of the plane at low viewing angles, 
they coincide at higher orientation angles (i.e., when seen 
as radio galaxy core emission). 

3.2. Radio Core Dominance 

In order to investigate the effect of viewing angle in 
populating the blazar envelope, it is necessary to have a 
way to measure or approximate the degree of misalign- 
ment in individual sources. The most readily available 
measure is the radio core dominance Rce; typical values 
range from —3 (for the most off-axis radio galaxies) to 
^ 3 (highly beamed quasars). As the extended emission 
is non-variable and unaffected by beaming, the depen- 
dence of Rce on the orientation angle 6 is directly due 
to the increased Doppler boosting factor at low angles. 
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Fig. 5. — Core dominance (Rce) versus synchrotron peak lu- 
minosity for BL Lacs, FSRQ, and radio galaxies in our sample. 
Color indicates the kinetic jet power (L]jin in ergs s~^), as calcu- 
lated from extended radio flux measurements, as in Figure |4] The 
dotted line has a slope equal to 1.2, the OLS bisector slope found 
for the whole sample. This slope is consistent with the expected 
relation between Lp^ak a-^d Lr due to beaming effects. 

The Doppler factor S is given by 
1 1 



(5 = 



r (1 - ^cosO) 



(3) 



where F is the Lorentz factor of the jet and /3 = v/c. The 
Doppler-boosted luminosity for any part of the jet SED 
is a function of S, 

L - LoSP (4) 

where L' is the unbeamed jet luminosity and p = 2 + a 
to 3 -f q; depending on the geometry of the jet model, 
with a equal to the e nergy spectral index at the fre- 
quency of interest (see lUrrv fc Padovanilll995() . Rce is 
thus a monotonically decreasing function of 9, though 
the exact shape is dependent on F, becoming steeper as 
increases. This does create a difficulty in comparing 
Rce values directly for sources which may have very dif- 
ferent bulk F values (such as, perhaps, very high and low 
power sources). However, if we assume that sources of 
the same jet power are similar in F, then we can com- 
pare at least these directly, and in general a source with 
a lower Rce value will be more misaligned. 

In order to investigate further how the sources of simi- 
lar jet power are related in FigureSJ it would be helpful to 
add the core dominance. We thus examine how the core 
dominance changes with peak luminosity and frequency 
separately in Figures [5] and [6l still continuing to bin by 
the jet kinetic power as before. This choice is partly be- 
cause it is practically difficult to add a 4th parameter 
in Figure m and also because, as suggested in Figure El 
the core dominance may have a different maximum (i.e., 
normalization), depending on the jet power, such that 
high-power sources have a lower maximum Rce than that 
of the lowest power bin, making it difficult to compare 



sources at different powers directly. This can be under- 
stood if we make the assumption that unbeamed radio 
power scales from the overall jet power as L^ = eLkin im- 
plying that jet radiative efficiency does not change much 
as the jet power increases. As we know Lkin=''^L(;^j (tak- 



ing 



0.6, see i )2.2p . this leads to the relation 



max Rce = (1 - l/l3) log Lki„ + p x log (2F) + ci (5) 

where we have let (5niax=2F, ci = log(eK^/''), and (1-1//3) 
~ —2/3. Assuming the last term is roughly constant, 
the maximum value of Rce is thus expected to decrease 
with increasing Lkin; even assuming reasonably different 
F values for high and low-power sources. For example, 
taking p=2, and F=15 for a jet of Lkin = 10^^ predicts 
a max Rce about 1.25 lower than for a jet with Lkin ~ 
10^2 and F=3. 

It is also expected that sources will separate by jet 
kinetic power in Figure [SJ as following from the above 
equations, we have for the expected relation between Rce 
and Lpcak: 



Rce 



aR 



log Lpcak - 



1 ^ o ) log Lki, 



-C2 (6) 



It is clear that binning by Lkin will indeed induce the 
'striping' in that figure through the second term, and 
that the slope in Figure [5] should be about 1.2 - 1.5 de- 
pending on the value of an (where we have assumed p = 
2 + a; see derivation in appendix IB|) . However, scatter 
could be induced by the remaining term, C2 = log(eK^' ") 
— (2+a_R)logL' j^j^/3 which is clearly determined by the 
relation between the unbeamed luminosity of a source 
at a given Lkin- If this quantity is not roughly constant 
with respect to Lkin, Figure [5] would be reduced to a 
multi-color scatter plot and the striping would be lost 
or degraded. It is also interesting to note that, like in 
Figure [21 there is a fairly uniform, sharp transition from 
blazars to radio galaxies in Figure [S] These arguments 
again suggest that beaming has a dominant effect on the 
observed (core) luminosity, rather than being due to a 
range of unbeamed core luminosities (L^, L' ^^j^) for a 
given Lkin- 

A mild decrease in synchrotron peak frequency with 
increasing angle to the line of sight is expected in the 
simplest case of a jet dominated by a single Lorentz fac- 
tor. The change in peak frequency as core dominance 
decreases is not significant for our sample of sources with 
Lkin greater than lO'^'*'^ ergs s~^, all of which have low 
synchrotron peak frequencies (these sources are therefore 
not shown in Figure JH]). For low-power sources (less than 
-[^q44.5 gj-gg g-i)^ however, there is a surprising pattern in 
the Rce — t'poak plane. As shown in Figure JH there ap- 
pears to be a clustering of sources with low peak frequen- 
cies and yet high core dominance, where Rce reaches 
up to ~3, and another location at high peak frequencies 
where Rce reaches ~ 1—1.5. Either locus could conceiv- 
ably be connected via a misalignment track to the radio 
galaxies at lower left, but the source distribution sug- 
gests that there is not a smooth continuum between LSP 
and HSP sources. We expect a similar monotonic rela- 
tionship between Rce and i^pcak as was derived above 
for the peak luminosity, but this is clearly not consistent 
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Fig. 6. — Core dominance (Rce) versus the synchrotron peak fre- 
quency for all sources with Ljjin below 10^^'^ (L300 below ~ 10^^) 
ergs s~^. While all low (kinetic) power radio galaxies have lower 
core dominance than blazars of similar power, there are apparently 
two locations of more aligned low-power blazars — one at low peak 
frequencies, with relatively higher Rce values (up to 3), and an- 
other at high peak frequencies, with slightly lower Rce values (~ 
1—1.5). While nearly all sources at high peaks (> 10^^ Hz) are 
BL Lacs (i.e., HBLs), the low-peaking aligned sources include a 
mixture of FSRQ and BL Lac spectral types. Color scale is the 
same as used in Figure [S] 

with what is shown in Figure [6l at least if sources at low 
power consist of a single population. 

In some cases the high core dominance values for the 
LSP sources at low powers could be due to systematic un- 
derestimation of L300 for these objects because of known 
problems with VLA measurements (incorrect modeling 
or use of an array configuration which might miss ex- 
tended flux). However this would require nearly every 
single source to be over-estimated by an order of magni- 
tude or more, and does not explain those LSP sources for 
which the steep spectrum is actually seen and matches 
the VLA estimate (such as BL Lac itself, noted in Figure 
li). 

4. DISCUSSION 

4.1. The Blazar Envelope 

iGhisellini et al.l ()1998D wer e able to reproduce the ob- 
served anti-correlations of the lFtassati et all (119981 ) blazar 
sequence using a simple one-zone leptonic jet model. 
They suggested that the increased radiative cooling of 
the electrons of high-power blazars lead to a particle en- 
ergy distribution with a break at lower energies (leading 
to the lower peak frequency). In support of the scenario 
of a continuous spectral sequence, our results in Figure S] 
clearly show a decreasing upper limit on Lpeak as fpoak in- 
creases. That sources near this limit are highly aligned, 
and those falling below are progressively misaligned is 
supported by the strong correlation between core dom- 
inance and peak luminosity shown in Figure [5] (which 
we note is much harder to appreciate without separating 
sources according to their radio galaxy or blazar classi- 
fication and by their jet kinetic power). However, this 



simple scenario does not provide an explanation for the 
rather noticeable 'L' shape to the source distribution in 
Figure |4] or the wide range of synchrotron peak frequen- 
cies for sources of similar low jet powers. While the for- 
mer might be the result of missing sources at intermedi- 
ate peak frequencies, the latter is a positive finding that 
cannot be explained away as a selection effect. 

4.1.1. Jet Structures in the Envelope 

IGhisellini fc Tavecchiol (|2008[ ) revisited the theoretical 
blazar sequence, linking the power of the jet and the SED 
to two fundamental parameters of the accretion disk — 
the black hole mass and the accretion rate. Following 
INaravan et al.l (jl997[) they assumed that sources with an 
accretion power (Lace) lower than a critical value Lcr ~ 
LEddfricr are radiatively inefficient accrctors with weak 
or absent broad line region, while sources with Lace > 
Lcr are characterized by efficiently radiating accretion 
disks with strong big blue bumps and powerful broad 
line emission (where LEdd is the Eddington luminosity 
and rhcr ^ 3 x 10"'^ — 10~^ a critical mass accretion rate 
in units of the Eddington accretion rate) . 

As discussed in ^3] we see in our results the manifesta- 
tion of two populations distinguished by typical jet struc- 
ture. The first population have non-decelerating jets 
characterized by a single Lorentz factor. These sources 
do not contain a high-peaking synchrotron component 
and the increase in peak frequency as the jet is aligned is 
modest, leading to a concentration of blazars with syn- 
chrotron peak frequencies ~ 10^^ Hz even when fully 
aligned. The sources in this population include all of 
the more powerful blazars (Lkin > 10"^^-^ ergs s~^) and 
many sources at lower powers. We include the latter be- 
cause of the high core dominance of many low-power, 
low-peaking blazars in Figure [5] The second population 
is only comprised of sources with Lkin below ^ lO^'*'^ 
ergs s~^ and have complex, decelerating jets. These 
lower-power sources exhibit a high frequency-peaking 
synchrotron component when aligned. 

The non-decelerating, low-peaking ('strong' jet) 
sources and decelerating, high-peaking ('weak' jet) 
sources form distinct populations which follow qualita- 
tively different de-beaming paths in Lpoak — t'poak space, 
exampled by track (A) and (B) in Figure |4l respectively. 
The two population scenario explains the appearance of 
an L-shaped upper boundary as a natural consequence 
of the vertical and horizontal paths by sources through 
the Lpeak — ^poak plane. It also explains the appearance 
of Figure [HI as low power sources may be either strong 
or weak, producing considerable spread in possible i^poak 
values as Rce increases. 

It is tempting to associate those weak sources (having 
Lkin < lO'''*-^ erg s~^) that we argue are characterized 
by velocity profiles (track B in Figure |4]) with sources 



havi ng Lg 



< Lc 



If we assume that Lt 



< Ln 



(e.g. iRawlings fc Saunders! Il991f ). this implies that the 
accretion power of these systems reaches at least up to 



10 



44.5 



erg 



which for rhr 



3 X 10- 



responds to a black hole mass of ~ 10 M (7), in agree- 
ment with black hole mass estimates (e.g., iBarth et al.l 
l2002f ). Similarly, sources with Lkin > 10*^'^ erg s~^, 
which are all of the powerful strong jet type (track A 
in Figure |4]), must have Lace > Lcr or equivalcntly 
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m > TTicr ^ 3 X 10^^, assuming that black hole masses 
do not significantly exceed ~ 10^ M0. 

In this framework, the low jet power sources (Lkin 
< 10^^'^ erg s~^) that yet seem to have a strong jet may 
be explained as sources with relatively smaller black hole 
masses, such that although their jets are low power, the 
mass accretion rate required to sustain them is m > rhcr- 
Consider a jet with Lkin = 10^^ erg s~^ and black hole 
mass 10* Mq. From the above, it is clear that the ac- 
cretion rate (in units of 10""^) m_3 = 0.7 Lkin, 44 ^g^ 
~ 7, exceeds the critical value placing this relatively low 
power jet in the family of 'strong' jets. 

4.1.2. Optical Spectral Type and the Envelope 

While th e ori ginal observational sequence of 
iFossati et al.l (|1998f ) implied a transition from quasars to 
line-less BL Lacs along a sequence from high to low bolo- 
metric luminosity, a firm distinction is not a necessary 
part of the theoretical scheme for a continuous blazar 
sequence. iGhisellini fc Tavecchiol ()2008D explain the 
possibility of blue FSRQ, formed when the dissipation 
radius for the jet exceeds the radius of the broad line 
region, as well as lower-power FSRQ formed when an 
efficient accretion forms around a relatively small black 
hole, (the explanation we proposed in ij4.1.1l is along the 
same lines). However, the traditional divide between 
BL Lacs and FSRQ, at an equivalent wi dth of 5 A may 



prove to be an essentially arbitrary one [G hiscUin i et al 
20101) . As discussed in IGeorganopoulos fc Marschei 
1998[ ). a line- less blazar may arise as a product of 
beamed emission out-shining the spectral lines, which 
would then emerge if the source were misaligned. A 
further complication arises when considering the spectral 
types of radio galaxies. High excitation radio galaxies 
(HERG) have prominent emission lines and presumably 
standard accretion disks accreting near the Eddington 
limit, while low excitation radio galaxies (LERG) have 
few or no lines and arc likely radiatively inefficient. 
While most FR II arc HERG and most FR I are LERG, 
there are known exceptions in both cases, includ ing the 
nearby HERG FR I Gen A (lEvans et al.l I2004D other 
possible 'FR I qua s ars' (iBlundeh fc Rawlingsl 120011: 
iHevwood et~all 120071: iChiaberge et al.l |2009[) . and an 
increasing pop ulation of LERG FR II at higher redshifts 

(up to z - O.S. lGarofalo et al.l[20Tol). 

An explanation is given by iGarofalo et al.l (|2010f ) for 
the presence of these mixed-type AGN based on black 
hole spin — namely, that central black holes with high 
retrograde spins and efficient accretion will produce pow- 
erful HERG FR lis, which slowly spin down to prograde- 
spinning, low efficiency FR Is. However, if the transition 
between spin state occurs significantly before or after the 
transition in accretion efficiency, mixed-state objects will 
occur. Interestingly, the spin-based unification scenario 
may explain the lack of an obvious continuum of blazar 
types, as in this scenario the transition objects (near zero 
spin) have much lower jet powers than rapidly spinning 
FR I or FR lis, thus in our typical luminosity-dependent 
plots there will be a 'missing population' of transition ob- 
jects, either showing up at much lower powers or missing 
due to selection effects. 

Ultimately, it is not clear what role the ionizing con- 
tinuum or line emitting regions has in influencing the jet 
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Fig. 7. — The ratio log (Lcore/Lx) versus synchrotron peak 
frequency for all 653 sources with well-measured synchrotron 
peaks (equivalent Orx is shown on the right scale, Orx = log 
(Lcorc/Lx)/7.68). This ratio is often used to estimate the location 
of the syn chrotron peak frequency when detailed SED information 
is lacking. ILandt &: BignaH 1(2003) use log (Lcorc/Lx) < 5 to indi- 
cate high synchrotron peak, however as shown, this ratio may only 
give a reliable upper limit to the frequency, and low log (Lcoro/Lx) 
values do not correspond exclusively to high synchrotron peaks. 



type. Optical spectral type is therefore not a part of our 
suggested strong/weak classification, which relies instead 
on the source morphology, core-dominance, and location 
in the i^poak— Lpcak plane. Thus, many of the low-peaking 
sources with high core dominance (i.e., strong jets) are 
BL Lacs, including BL Lac itself and all but one of the 
nine hot-spot BL Lacs found in the MOJAVE sample 
(jKharb et al.ll2010f ). Similarly, several quasars appear at 
higher peak frequencies (i.e., as weak jets) in the UEX 
sample. 

Lastly, we comment on the possible presence of high- 
frequency-pea king FSRQ (HFSRQ) as claimed by sev- 
eral authors (ILandt et al.l 120061: iPadovani e\~^ 120071: 
ILandt et al.ll2008t ). While the presence of quasars with 
low bolometric luminosities and high peaks (essentially 
having SEDs typical of HSP BL lacs) are consistent with 
expectations based on the up dated theoretical sequence 
(jGhisellini fc Tavecchiol l2008| ). many of the objects dis- 
cussed do have higher jet kinetic powers typical of FR 
lis. However, using low core radio to X-ray flux (log 
Lcoro/log Lx < 5) or high optical to X-ray ratio {aox) 
is not a reliable way to find sources with high t'pcak- To 
illustrate this, in Figure [7] we plot log (Lcorc/Lx); mea- 
sured at 5 GHz and 1 keV, against the known peak from 
our SED fits (based on much more extensive data than 
a two-point spectral index) for the 653 sources in our 
TEX+UEX sample. While a few quasars exist in the 
HSP range > 10^^ Hz, these sources appear otherwise 
very similar to HSP BL Lacs in their SEDs. Further, 
many sources below the log (Lcore/Lx) < 5 limit clearly 
have low peak values, often with a hard X-ray spectrum 
indicating a non-synchrotron origin. So far no unambigu- 
ous evidence of any blazar with a high peak frequency 
and high luminosity (above the typical HSP zone) has 
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TABLE 4 
Overview of Strong and Weak Jet Classes 



Strong 



Weak 



Definition 

Optical Type 
Synchrotron Peak 
Typical L300 
Typical L^in 
Radio Morphology 
Accretion Type 



More powerful, single-component jet (single char- Lower-power jet with faster deceleration, multiple 
acteristic F), remains relativistic on large scales emitting regions (multiple F) 



mostly FSRQ 

All high-power LSP, some low-power LSP 

lO-lO - 10"'5 ergs s-i 

10-13.5 _ ;^o46 ergs s-i 



mostly BL Lacs 

All ISP, HSP, some low-power LSP 

1038 _ 10*2 ergs s-i 

10*2 - 1044-5 ergs s^i 



FR II; CoUimated jets, hot spots, edge- brightened FR I; Less coUimated, edge-darkened, no hot spots 
Efficient, likely standard disk Inefficient, highly sub-Eddington 



been found fsee lMaraschi et al.|[2008HLandt et al.ll2008D . 

4.1.3. Intermediate and Low-peaking Blazars 

According to the proposed two-population scenario, all 
ISP and HSP sources belong to the weak jet parent sam- 
ple along with FR I radio galaxies. As blazars, they move 
from a population of HSP to ISP and finally LSP sources 
as they are misaligned. Contrary to the int erpretation 
of a continuous sequence (jFossati et al.|[l998l ). then, ISP 
blazars should have lower core dominance, lower overall 
luminosities, lower Doppler factors, and yet similar Lkin 
when compared with HSP sources. However, it is difficult 
to rule out a missing population of intermediate sources 
which may fill in the gap between low and high peaks in 
Figure [6] or in the center of Figure ID These less popu- 
lated regions are not due to the fiux-limited selections of 
the original master sample (see below) . If these sources 
are in the master sample, the SED sampling would have 
to be much lower, and for such relatively bright objects, 
this seems unlikely. 

Weak jet sources at even greater angles to the line of 
sight should appear as LSP blazars, though they may not 
be recognizable as such (but rather as radio galaxies). 
Regardless, these should be distinguishable from strong 
jet LSP sources based on core dominance, morphology, 
extended radio power, and polarization signatures. How- 
ever without detailed VLBI and VLA studies to measure 
these attributes, confusion among the two LSP popula- 
tions is likely to occur. From our hypothesis that Rce 
must decrease with peak frequency for weak-jet sources, 
we expect most of the low-power LSP blazars shown in 
Figure H] to belong to the parent sample of strong jet 
sources, along with FR II radio galaxies. This connec- 
tion is suggested by the equal or higher core dominance of 
most of these sources compared with the HSP blazars in 
the same Lkin range. It is also supported by the detailed 
observations we do have from VLBI studies of some of 
our objects, such as the nine LSP BL Lacs in the MO- 
JAVE sample with hotspots similar to those typically 
seen in high-powered quasars and associated with FR II 
morphology (JKharb et al.ii20To[ ). The MO J AVE blazars 
(almost entirely LSP sources) were also found to have an 
intrinsic, unbeamed luminosity function consistent with 
that of FR II radio galaxies, whet her BL Lacs were in- 
cluded or not (jCarafc ListeHl2008f ). 

We summarize the strong and weak jet characteristics 
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Fig. 8. — The blazar envelope for the UEX sample of 437 sources 
with only upper limits on Lkin. Upper limits were found by extrap- 
olating from the lowest frequency observation of core flux using a 
spectral index a = 1.2. The same color scale is used as in Figure 
[3] We clearly see a similar 'L' shape as is suggested in the flgure 
for TEX sources. Also shown are the three zones discussed in i|4.2l 

as suggested by the current work in Table 31 

4.2. Examining the Sample Bias 

Our effort is concerned with those sources which can be 
placed on the Vpcak ~ Lpcak plane with confidence, which 
is primarily a function of SED sampling. It is important 
to examine both the selection of the entire sample and 
the nature of the discarded objects for evidence of any 
bias which might produce the current findings. 

One possible selection bias is a 'windowing' effect due 
to the low frequency radio measurements required to iso- 
late the core and extended components. On one side, 
very core-dominated sources are likely to be missed when 
using only the spectral decomposition method, as these 
might not exhibit a steep spectrum except below very low 
frequencies (10 — 100 MHz) which arc not frequently ob- 
served (and thus fail the two-component test). On the 
other side, a very de-beamed source may not qualify as 
a flat-spectrum source if the Across is much above typ- 
ical sample definition frequencies of a few GHz. The 
combined effect could in principle lead to the selection 
of sources for which the core and extended fluxes are 
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by constraint rather similar. The inclusion of VLA map 
data, which allows us to find extended flux densities even 
in very core-dominated sources, should mitigate the first 
issue. As to the second, as we show in Figure [21 most 
blazars arc core dominated above the crucial window of a 
few GHz, and we detect a few above this limit. Between 
the blazars and radio galaxies in this figure there is no 
large gap indicating a missing population, and we do not 
find any reasons to suspect that the windowing effect is 
a problem for the results of this paper. 

4.2.1. Upper Limits 

The over 400 sources in the UEX sample are shown 
in Figure [H with the color corresponding to the upper 
limit of the jet kinetic power. The overall shape of the 
envelope appears qualitatively similar to Figure |4l 

We tested areas of the fpcak~Lpcak plane using simula- 
tions to check how the sensitivity to sources in the plane 
is affected by our overall combined sample selection (i.e., 
the combined selection of the samples in Table [1]) . We 
moved ~100 template SEDs drawn from our database to 
various peak locations in each zone shown in Figure |8] 
and simulated cuts using the flux limits of our individual 
samples. In zones 1 and 2, we did not flnd any difficulty 
detecting the simulated sources out to a reasonable red- 
shift (z ~2). In zone 3, however, our sample would likely 
only include sources out to z ~ 0.2, which is not much less 
than the typical redshift of the HSP sources just above 
this region (z ^ 0.25). Wc conclude that the particular 
lack of ISP sources with Lpoak intermediate between LSP 
and HSP sources is not due to the initial selection of our 
sample. For at least zones 1 and 2, which would contain 
very bright sources, it is unlikely that a large number 
of sources would exist there, but lack the SED sampling 
necessary to make the more subjective cuts leading to 
our final sample. 

4.2.2. Remaining Sample 

We now examine the sources which wc were unable to 
include due to inadequate SED sampling (typically fewer 
than ~ 10 SED data points). Examination of the colors 
(radio-optical and optical-X-ray indices) for these sources 
indicates that they have values consistent with the well- 
characterized sample of ^ 600 TEX-I-UEX sources. Un- 
surprisingly, those with X-ray detections mostly popu- 
lated the 'HSP zone' of the a^o — chox plane, indicating 
that a large population of the excluded sources are likely 
higher redshift HSPs and ISPs. 

The optical and X-ray luminosity distributions for the 
excluded sample as a whole are indistinguishable (by KS 
test) from the TEX-fUEX sample. In all three bands 
(and in particular the radio), those sources detected in 
the X-rays tended to have much lower luminosities, and 
are therefore unlikely to populate the so-called 'forbid- 
den zone' of high synchrotron peak luminosity and fre- 
quency. Those not detected in X-rays had a distribution 
of aro typical of LSP sources, but actually tended to have 
higher luminosities than our TEX-fUEX sample. How- 
ever, non-detection in X-rays and higher redshifts in gen- 
eral probably explains their absence in the TEX+UEX 
sample. Overall, none of the sources excluded from the 
TEX-I-UEX sub-sample appear to be substantially dif- 
ferent from those in it; however, only additional multi- 
frequency observations of these sources can confirm this. 



5. SUMMARY AND CONCLUSIONS 

In this work we have analyzed the multi-frequency 
spectrum of several hundred blazars in order to accu- 
rately measure the synchrotron peak frequency and lu- 
minosity, and performed spectral analysis of the low- 
frequency emission in order to estimate the intrinsic jet 
power by a scaling of the 300 MHz extended radio power. 
The connection of blazars and radio galaxies in the same 
plane amounts to an extension of the original blazar 
sequence to a blazar envelope, in which highly aligned 
blazars form an upper limit in the :^pcak — Lpcak plane, 
while radio galaxies populate the lower left. We flnd that 
the jet power does not absolutely determine the i/pcak 
of a maximally-beamed source. Instead, there seems to 
a break-down at low kinetic jet power, such that these 
jets are capable of having either low or high synchrotron 
peaks when aligned with the line of sight. 

Our results suggest a shift away from the continuous 
blazar sequence interpretation, but confirmation waits 
for ever more data. We summarize the most important 
points here: 

• Among our sample of over 600 blazars with well- 
sampled blazar SEDs, none with high peak frequen- 
cies and high luminosities are found. This void may 
be explained by the orig inal arguments advanced 
by iGhisellini et al.l ()1998D , that higher external ra- 
diation fields in high-power sources lead to greater 
cooling and thus lower break/peak energies. 

• However, it is not clear that there is a continuum 
of SED types with jet power. The source power 
alone does not appear to determine the placement 
of a source along the aligned sequence. Instead, we 
suggest that radio-loud AGN can be divided into 
two populations: low-synchrotron-pcaking sources 
without significant velocity gradients or deceler- 
ation (strong jets) and high-synchrotron-peaking 
sources which do have such structure (weak jets). 
The attributes of these samples are summarized in 
Table 12 

• It appears that below a certain jet power, 
sources can either have strong or weak jets, sug- 
gesting another parameter besides kinetic power 
is responsible for the jet structure. As in 
IGhisellini fc Tavecchiol (|2008[ ). we suggest that this 
parameter is the accretion rate m in Eddington 
units with the two families separated at m ~ 
rhcr ^ 3 X 10"'^ — 10~^, with sources having 
TO < rhcr exhibiting rad iatively inefflcient accre- 



tion (jNaravan et al 



L g rad i 

MM) 



• A key prediction of the two-population scenario is 
that ISP sources primarily consist of misaligned 
HSP objects, as suggested by their slightly lower 
Lpcak values and similar Lkin to HSPs. Observa- 
tions of ISP sources with synchrotron peak values 
and core dominance values between those typical 
of the LSP and HSP sources (currently missing) 
would support a continuous sequence interpreta- 
tion instead. 

• Observations of highly core-dominated, low- 
peaking BL Lacs (such as BL Lac itself), including 
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several BL Lacs exhibiting FR Il-likc morphology, 
leads us to speculate that these sources belong to 
the strong-jet parent population, along with FR 
II radio galaxies. Likewise, some ISP and HSP 
sources with quasar-like spectra in our sample be- 
long to the weak population. Thus the spectral 
types associated with strong and weak jets are 
mixed. 



The low peak frequencies seen in our FR I radio 
galaxies supports the interpretation of a horizon- 
tal movement in the i^peak — Lpoak plane for weak 
sources which supports a model with velocity pro- 
files in the jet and leads to our two-population 
scheme. Observations of high synchrotron peaks 
in radio galaxies would be very difficult to explain 
under this scenario. 



The authors wish to thank the anonymous ref- 
eree for helpful comments, and Katherine Blundell 
for a discussion which lead to the improvement of 
this manuscript. This research has made use of 
the SIMBAD database, operated at CDS, Stras- 
bourg, France (http://simbad.u-strasbg.fr/simbad/), 
and the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Ad- 
ministration (http://nc dwww.ipac.calte c h.edu/ ). We 
utilized package np (Hayfield & Racine 2008|) in the 
R lang uage and environment fo r stat istical com- 
puting ()R Development Core TeamI |2010[ ). GF ac- 
knowledges support from NASA grants NNG05GJ10G, 
NNX06AE92G, and NNX09AR04G, as weh as SAO 
grants G03-4147X and G05-6115X. MG acknowledges 
support from the NASA ATFP grant NNX08AG77G and 
NASA FERMI grant NNH08ZDA001N. The MOJAVE 
project is supported under National Science Foundation 
grant 0807860- AST. 



APPENDIX 



A. PHENOMENOLOGICAL BLAZAR SED MODEL 

The model was adapted from iFossati et all (J1997I ) and used to fit the average blazar SED with seven parameters, 
six of which were free to vary and which were optimized using an SA routine. 
Parameters: 



T E 



•^peak 

/3 = 



32-46 

1000 

1-4 

3-9 

-0.4-0.8 

10-20 

0.6 - 2.0 



To build the SED, the following equations were combined, where x = log(;7). Vector quantities are noted with arrows 
for clarity, and luminosities are in log. 



2x; 



(Al) 



-'peak 



li (Xpcak -X;- 9.698) + Lr + \-li 



(A2) 



^synch 




^damp 



when X < Xpeak 
else 

when X < (xpcak - Xj) 
else 



(A3) 



(A4) 



^ic 



K,.x + 7.685 + (1 - a.^) (x - 17.383) 



(A5) 



Finally, we have the complete SED: 



L = log 10^=^"<='' + 10^i<^ 



(A6) 
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B. RADIO CORE DOMINANCE 
Equation |6] can be derived starting from the definition of log core dominance ratio Rce: 

Lr. 

-Loxt , 



RCE = log 



(Bl) 



where Loxt and Lr are the luminosities of the extended and core components taken at the same frequency, and the 
beaming equations for the luminosity for the radio and synchrotron peak frequency are simply: 



T , — T ' A^'P"""' 
iJpcak — ^pcak"poak 



(B2) 
(B3) 



where 5r and ^poak are the Doppler factors for the jet regions dominating the emission at those frequencies, and the 
exponents take values p = 2+a or 3+a, depending on the jet model, and with a the spectral indices in the band of 
interest. By taking advantage of the relationship Lkin = i-^^cxt (see ij3.2p we can rewrite it as: 



RcE = log 






(B4) 



If we assume that Lp = eLkin (see ^ ^3.2(1, a nd that the bulk Lorentz factors T for the radio and the peak emission are 
the same, such that S^ ~ i5pcak, using lB3l wc obtain: 



RcE = log eLk 



Letting pa = 2 + aR and ppoak = 2 



^pcak 



R-CE — 



an. 



I-'pcak 
peak 



log 



= 3, we have equation [51 



log L 



peak 



C2 = log(e«:i/^) 



CtB, 



Lkir 
K 



1//3 



1 - - ) log Lkin + C2 



logL; 



peak 



(B5) 

(B6) 

(B7) 
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